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Abstract 

A measurement of the spin alignment of charged D* mesons produced in 
continuum e + e~ — ► cc events at = 10.5 GeV is presented. This study 
using 4.72 fb _1 of CLEO II data shows that there is little evidence of any D* 
spin alignment. 
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I. INTRODUCTION 



There have been numerous theoretical JT|-[7| and experimental [PHT5H studies of the frag- 
mentation of heavy quarks. The energy distribution and flavor dependence of heavy quark 
hadronization have been modeled by fragmentation functions. The role that spin plays in 
the hadronization process is still being investigated and is not well understood at this time 
|T6^f22| . To increase the understanding of this role, a precise measurement of the probabilities 
of a meson being directly produced in each of the available spin states is needed. 

At CLEO, the fragmentation of charm quarks can be analyzed by making measurements 
of primary hadrons containing charm quarks from continuum e + e~ annihilations. CLEO has 
previously published results of charmed meson energy distributions || as well as the spin 
alignment of charged D* mesons |18|]. In this paper an updated measurement of the charged 
D* spin alignment using the entire CLEO II dataset is presented. 



II. POLARIZATION, ALIGNMENT, AND P v 

According to the quark model, a meson is composed of two spin | valence quarks that can 
combine to form four spin states in the absence of orbital angular momentum, i.e. four S- 
wave states. Writing these in the basis of total angular momentum, J, and its z-component, 
J z , they are the vector states |1,1 >, |1,0 >, |1, — 1 >, and the pseudoscalar state 1 , >, 
where the z-direction can be arbitrarily chosen. The probability of an S-wave meson being 
produced in a vector state is often described by the ratio Py defined as 

Pr v 



V + P 

where P and V are the respective probabilities of the meson being created in the pseudoscalar 
and vector states. 

The helicity formalism is useful in the context of describing the angular distributions and 
correlations in the production and decay of particles with non-zero spin. For a particle with 
momentum p, the helicity is defined as 

which in the case of a spin-1 particle is just the z-component of the spin when the z-direction 
has been chosen as the flight direction of the meson. The helicity density matrix is often 
used to organize information about the spin of a particle. The diagonal elements of this 
matrix p\\, with J2\Pxx — 1> represent the probability that the particle has helicity A. 

Simple statistical expectations are that all helicity states of a spin J particle are equally 
populated, but production and fragmentation dynamics can lead to either polarized or 
aligned particles. A system of particles is polarized if there is a net angular momentum, 
i.e. px\ 7^ P-a-a for some helicity A, and it is aligned if there is a nonuniform population of 
states, but paa = P-a-a f° r an A. Since the production and fragmentation processes in this 
analysis conserve parity and the CESR beams are unpolarized, it is expected that the D* 
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mesons from e + e~ —>■ 7* — > cc are unpolarized, but it is possible for the D* mesons to be 
aligned. 

To measure the spin alignment of a vector meson, the angular distribution of its decay 
products is analyzed, but because the angular distributions of the A = 1 and A = — 1 states 
are degenerate, the values of pn and p_i_i cannot be distinguished and only one variable, 
e.g. poo — 1 — P11 — P-i-i) is accessible. From the definition above, the vector meson is 
aligned if poo differs from 1/3. For the case of a vector meson decaying to two pseudoscalar 
mesons, the angular distribution can be written 

W(cos9) = ^[(1 - poo) + (3p 00 - 1) cos 2 #] (3) 

where 6 is defined as the angle of a daughter pseudoscalar in the parent vector meson rest 
frame, with respect to the direction of motion of the parent vector meson in the rest frame 
of the production process. In our case, the production rest frame of a D* directly produced 
in charm fragmentation from e + e~ annihilation coincides with the lab frame. 
By using the variable 

3poo — 1 /a 

« = 1 , ( 4 ) 

1 - Poo 

the angular distribution can be expressed as 

W(cos6) = N(l + a cos 2 6) (5) 

where N is a normalization factor equal to 3/ (6 + 2a). The value of a can range between —1 
and +00, where the angular distribution would be isotropic if a = 0, proportional to sin 2 6 
if a = — 1 and proportional to cos 2 9 if a = 00. 

Whereas the naive statistical expectation is that all four 5-wave meson states are created 
in equal proportions, i.e. p\\ — | (a — 0) and Py = 0.75, there are other models that have 
been presented where the alignment and Py vary as a function of momentum P3"|Z1|. Heavy 



quark symmetry predicts that vector mesons containing a single heavy quark are produced 
unaligned, but there have been suggestions that the value of Py may depend upon the mass 
difference of the vector and pseudoscalar mesons P,p5f . It has also been suggested that Py 
is directly related to the spin alignment [^] and in the previous CLEO D* spin alignment 



analysis, a value for Py was calculated using this relationship [18]. However, the validity of 
the statistical model is assumed when deriving this relationship. We feel that a determination 
of Py for D* mesons warrants an independent measurement which is the topic of a current 
CLEO analysis. The results of the Py analysis will be presented in a future paper. 



III. DETECTOR AND EVENT SELECTION 

The CLEO II detector is a general purpose charged and neutral particle detector and is 



described in detail elsewhere [p7 |. The dataset used in this analysis consists of 3.11 fb _1 of 
data collected at the T(4S) resonance and 1.61 fb" 1 of data collected about 60 MeV below 
the resonance. This corresponds to approximately 5 million continuum cc events. 

The D* + in this analysis is required to decay through the channel D* + — > D°tt + with the 
D° decaying either through the mode D° —* K~n + or D° — > K~7t + tt° (inclusion of charge 
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conjugate modes is implied throughout this paper). The 7r + in the D* + decay is kinematically 
limited to having a momentum less than 456 MeV/c in the lab frame of reference, and is 
referred to as the "slow" pion. 

All tracks used in this analysis are required to have an impact parameter within 5 mm 
of the interaction point in the plane transverse to the beam pipe and within 50 mm in the 
direction of the beam pipe. Tracks are also required to have a momentum less than 6 GeV/c 
and an r.m.s. residual less than 1 mm for their hits. Particle ID is not used since there is 
little gain for this particular analysis and it introduces the possibility of additional systematic 
errors. For a pair of photons to be considered as a candidate 7i~°, they must each have an 
energy of at least 100 MeV, be within the "good" barrel of the detector (| cos Odetector I < 0.71), 
have a shower shape in the crystal calorimeters consistent with that of a photon, combine 
to be within 20 MeV/c 2 of the neutral pion mass, and have | cos^l < 0.9, where # 7 is the 
decay angle of the photon in the ir° rest frame, with respect to the 7r° direction of motion in 
the lab frame. 

For the D° — > K~tt + mode, the D° is reconstructed by taking all possible pairs of 
oppositely charged tracks in an event, assigning the kaon mass to one and the pion mass 
to the other (or vice versa), adding their four-momenta, and then calculating the invariant 
mass. The D* + is reconstructed by adding the four-momentum of a candidate slow ir + in 
the event to the four-momentum of the candidate D°. The mass difference, AM, between 
the candidate D° and D* + is required to be within 2.5 MeV/c 2 of the world-average mass 
difference of 145.42 MeV/c 2 @. 

The D° is spinless and the decay products have an isotropic angular distribution. How- 
ever, due to the jet-like nature of continuum events, the background from random track 
combinations tends to have cos0x ~ —1, where 4>k is the decay angle of the K~ in the D° 
rest frame, relative to the D° motion in the lab frame. A requirement that cos0^; > —0.9 is 
added to improve the signal-to-background ratio. 

For the D° — > K~h + 'k mode, the four-momentum of a candidate 7r° is added to the 
four- momenta of two oppositely charged tracks to form candidate D 0, s in the event. Mass 
difference and kaon decay angle requirements are the same as described above. 



IV. FITTING 

To test models that predict that the alignment varies as a function of the momentum of 
the D* + , the data are broken up into six x + bins in the range 0.25 to 1.0, where 
Lorentz-invariant variable defined as 

x+ = P(D*) + E(D*) 

P (D*) + E (D*)' K ' 

1 maxy^ ) ~ J -'max\ J -' ) 

where E max = E beam , P max = E% eam - M% t+ and M D *+ is the world- average value for the 
mass of a D* + . 

For each x + range, a sideband subtraction is performed. The sideband region is from 
9 MeV/c 2 to 12 MeV/c 2 above the mean of the AM peak and the ratio for the sideband 
subtraction is determined by fitting the data with a bifurcated double Gaussian for the signal 
plus a background function A+ B (AM) 1 / 2 +C (AM) 3 / 2 and integrating the background shape 
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for the signal and sideband regions. The fits used to determine the sideband ratios are shown 
in Figures [l] and |[ 

The sideband-subtracted M(Ktt) data is fit for each x + bin with a double Gaussian for 
the signal region plus a first-order polynomial background. [] Each of these x + bins is broken 
up into five equal cos 9 bins, where 9 is the angle defined in Section II. To prevent the fitted 
signal shape from having large fluctuations due to lower statistics in some cos 9 bins, the D° 
mass peak is fit for each cos 9 bin with the ratios of areas and widths of the double Gaussian 
fixed to those found when fitting the mass peak in that momentum range for the entire cos 9 
spectrum. 



V. EFFICIENCIES 

It is important to understand the relative efficiencies of detecting a D* + in the various 
cos 9 bins. In the lowest momentum bins, for example, the efficiency decreases as cos# 
approaches one because of the increased difficulty in measuring the track of a slow pion that 
is emitted in the direction opposite the D* direction in the lab frame. Detection efficiency as a 
function of x + and cos 9 was measured by analyzing Monte Carlo data with a GEANT-based 
detector simulation. 

Monte Carlo events were generated using the Lund Jetset 7.3 program, where the e + e~ 
annihilation was required to result in a cc pair with one of the charm quarks hadronizing to 
a D* + that decays to D°n + with D° — > K~7t + (tt°), while no constraints were placed on the 
other charm quark. The D* mesons were produced such that their decay to D°7r + had an 
isotropic angular distribution in the rest frame of the D* + . 



VI. RESULTS 

The fits of the sideband subtracted M(Ktt) and M(Knii ) distributions for all scaled 
momentum ranges are shown in Figures [| and |J0 The efficiency-corrected angular distribu- 
tions for both decay modes were combined in each x + bin with a weighted average and are 
shown in Figure [5], where they have each been normalized to unit area and fit with Eq. (5). 

The values of a resulting from these fits as well as the fits for each of the two decay modes 
treated separately are listed in Table |. Figure |5| shows the combined results for a plotted 
as a function of momentum as well as the theoretical curves suggested by Suzuki |S| and 



Cheung and Yuan \24\. Table || lists the values of poo as calculated from the measurement 



of a for each scaled momentum bin. Averaging the cos 9 distributions over all momenta and 
then fitting gives a value a = —0.028 ± 0.026, corresponding to p 0Q = 0.327 ± 0.006. 



The highest x + bin is fit with a second order polynomial for the background since the background 
is not well represented by a straight line. 

2 Only the highest five momentum bins were used for the D° — > Kttt: mode due to the small 
number of signal events and low signal-to-noise ratio in the lowest x + range. 
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3110198-004 




160 



AM = M (Ktttt)- M (Ktt) (MeV / c ) 



FIG. 1. D* — D mass difference for the D° — ► i^7r decay mode for the six x + ranges a) 
0.25 < x+ < 0.45, b) 0.45 < x+ < 0.55, c) 0.55 < x+ < 0.65, d) 0.65 < x+ < 0.75, e) 
0.75 < x + < 0.85, f) 0.85 < x + < 1.0. The solid squares are the data points and the solid 
line is the fitting function as described in Section IV. The hatched area is the signal region while 
the cross-hatched region is the sideband. 
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FIG. 2. D* — D mass difference for the D° — > Ktttt decay mode for the five x + ranges 
a) 0.45 < x+ < 0.55, b) 0.55 < x+ < 0.65, c) 0.65 < x+ < 0.75, d) 0.75 < x+ < 0.85, e) 
0.85 < x + < 1.0. The solid squares are the data points and the solid line is the fitting function as 
described in Section IV. The hatched area is the signal region while the cross-hatched region is the 
sideband. 
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FIG. 3. M(Ktt) after sideband subtraction for the six x + ranges a) 0.25 < x + < 0.45, 
b) 0.45 < x+ < 0.55, c) 0.55 < x + < 0.65, d) 0.65 < x+ < 0.75, e) 0.75 < x + < 0.85, f) 
0.85 < x + < 1.0. The solid squares are the data points and the solid line is the fitting function as 
described in Section IV. 
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FIG. 4. M(K7nr°) after sideband subtraction for the five x + ranges a) 0.45 < x + < 0.55, b) 
0.55 <x + < 0.65, c) 0.65 < x + < 0.75, d) 0.75 < x + < 0.85, e) 0.85 < x + < 1.0. The solid squares 
are the data points and the solid line is the fitting function as described in Section IV. 
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TABLE I. Values of a for different momentum ranges. First error given is statistical, second 
is systematic. The last column is the confidence level of the fit for the combined values of a. 
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FIG. 5. Normalized cos 9 distributions in the six x + ranges for the D° -> R--k + and 
D° -► K- 7T + 7r° decay modes combined. The solid squares are the efficiency-corrected yields 
for each cos# bin in the specified x + range. These distributions are fit with the function 
W(cos9) = 0.4AT(1 + a cos 2 9), where the factor of 0.4 is the bin width and N = 3/(6 + 2a). 
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0.25-0.25 0.40±0.07±0.07 

0.45-0.55 0.31±0.03±0.01 

0.55-0.65 0.33±0.01±0.01 

0.65-0.75 0.30±0.01±0.01 

0.75-0.85 0.35±0.01±0.01 

0.85-1.0 0.35±0.02±0.01 

TABLE II. Values of poo for different momentum ranges. The first error given is statistical, 
the second is systematic. 
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FIG. 6. The values of a for each momentum bin are represented by the solid squares. Errors 
shown are the statistical and systematic errors added in quadrature. The solid line represents the 
statistical model. The dotted line represents the function predicted by Suzuki [22]. The dashed line 
is the function of Cheung and Yuan [23]. 
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Collaboration \fs (GeV) a poo 

HRS 29 0.18 ± 0.08 0.371 ± 0.016 

TPC 29 -0.14 ± 0.17 ± 0.03 0.301 ± 0.042 ± 0.007 

SLD 91 0.019 ± 0.378 ± 0.582 0.34 ± 0.08 ± 0.13 

OPAL 91 0.33 ± 0.11 0.40 ± 0.02 

CLEO 1.5 10.5 0.08 ± 0.07 ± 0.04 0.351 ± 0.015 ± 0.008 

CLEO II 10.5 -0.028 ± 0.026 0.327 ± 0.006 

TABLE III. Results for a and poo found by various collaborations. 

Similar analyses have been done by the HRS, TPC, SLD and OPAL collaborations 
j^ . [l7|^T| , |20|] , as well as by CLEO using a previous data set fllSfl . The average values of 



a and poo m each study are presented in Table [TTT]. 



VII. SYSTEMATIC ERROR 

Many possible sources of absolute systematic uncertainty, such as the overall track-finding 
efficiency, do not have a significant effect on this analysis because the extraction of a in each 
momentum range involves only the relative comparisons of the same measured quantity, 
namely the yield of the D° decays, in the different bins of cos 9. The remaining sources of 
uncertainty will therefore be related to extracting the yield and the efficiency as a function 
of cos 9. The effects of the various sources of systematic error are shown in Figure |7| while 
the methods used to determine these errors are described below. 

The Monte Carlo contribution to the systematic error was accounted for by including 
the error in the Monte Carlo efficiencies in the calculations of the yields. To investigate the 
systematic error associated with the fitting function, the analysis was done using a single 
Gaussian rather than a double Gaussian to fit the signal peaks. Likewise, to investigate the 
systematic error associated with the choice of range for the sideband subtraction, the analysis 
was done using a sideband region from 6 MeV/c 2 to 9 MeV/c 2 above the nominal D* — D 
mass difference rather than 9 MeV/c 2 to 12 MeV/c 2 above the nominal value. The effect of 
the mass difference requirement was investigated by constraining the mass difference to be 
within 1.25 MeV/c 2 of the PDG value rather than 2.5 MeV/c 2 . The systematic effects of 
the cos 9 binning were studied by using six equal cos# bins rather than five. The differences 
between the resulting values of a and the central value were all summed in quadrature as 
an estimate of the systematic error and are included in the error bars shown in Figure |6|. 

A small linear component in the angular distribution can easily be seen in Figure [5] for 
the range 0.65 < x + < 0.75. This is most likely due to a slight inaccuracy in the efficiency 
correction from the Monte Carlo data. The data in Figure |5| were fit with a straight line 
added to Eq. (4) as a check and the difference in the fitted values of a was negligible. 



VIII. INTERPRETATION OF RESULTS 

We have measured the spin alignment of all D* mesons produced in e + e~ — > qq interac- 
tions at yfs = 10.5 GeV. Although the details of the analysis ensure that the measured D* 
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FIG. 7. The results from the systematic error studies for the six x + bins a) 0.25 < x + < 0.45, 
b) 0.45 < x+ < 0.55, c) 0.55 < x+ < 0.65, d) 0.65 < x+ < 0.75, e) 0.75 < x+ < 0.85, f) 
0.85 < x + < 1.0. The different symbols represent the resulting values of a using the modifications 
in the analysis procedure as described in Section VII. 
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does not come from a decaying B meson, we cannot determine any other details about the 
production hierarchy. From a theoretical standpoint, we are particularly interested in the 
D* mesons that are produced directly in the e + e~ collision, but we cannot distinguish these 



from secondary £)*'s resulting from decays of charm mesons with L > 0. [29 31 



The most prominent excited charm mesons, which are commonly referred to as D** 
mesons, consist of a charm quark and a light anti-quark with relative orbital angular mo- 
mentum L = 1. They are categorized into four states with spin-parity J p = + , 1 + , 1 + , 
and 2 + . A + state decay to D*n is forbidden due to spin-parity conservation while other 
D* modes are expected to be suppressed. When a 2 + state decays through a D* channel, it 
can only produce a D* meson with a helicity of ±1 in the 2 + rest frame, while the 1 + states 
only decay through D* channels and favor a helicity of in the 1 + rest frame. From the 
measurements available [52"|,33|, we estimate that 16-20% of D* mesons observed at CLEO 
could be daughters of a D** meson, not including the contribution from D** mesons. 

Although the favored helicities of _D*'s from the decays of 2 + and 1 + charm states partially 
cancel, it is probable that these _D*'s are aligned in their production rest frame, i.e. the rest 
frame of the parent D**. It is expected that any effect would be most noticeable for the 
highest x + bins which has the largest correlation between the D* 4-momentum in the lab 
frame and the D* 4-momentum in the D** rest frame. If the 4-momenta in the two reference 
frames are uncorrelated, as tends to be the case for the lower x + bins, any alignment of _D*'s 
from _D**'s would not be noticeable in the lab frame. 

Due to the current lack of information about the production and decay of P-wave charm 
meson states, we can only state that D** decays could have a significant effect on this D* 
spin alignment measurement in at least some of the x + bins. 



IX. CONCLUSION 

This analysis is the most precise measurement of the spin alignment of D* + mesons to 
date. The data, without any corrections for D** effects on the measurements, agree well 
with the statistical model expectation that the J z = state has a | probability of being 
populated. 
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